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We formulate and investigate the simplest version of time-optimal quantum computation theory 
(t-QCT), where the computation time is defined by the physical one and the Hamiltonian contains 
only one- and two-qubit interactions. This version of i-QCT is also considered as optimality by 
sub-Riemannian geodesic length. The work has two aims: one is to develop a t-QCT itself based 
on physically natural concept of time, and the other is to pursue the possibility of using t-QCT 
as a tool to estimate the complexity in conventional gate-optimal quantum computation theory 
(p-QCT). In particular, we investigate to what extent is true the statement: time complexity is 
polynomial in the number of qubits if and only if so is gate complexity. In the analysis, we relate 
t-QCT and optimal control theory (OCT) through fidelity-optimal computation theory (/-QCT); 
/-QCT is equivalent to t-QCT in the limit of unit optimal fidelity, while it is formally similar to 
OCT. We then develop an efficient numerical scheme for /-QCT by modifying Krotov's method 
in OCT, which has monotonic convergence property. We implemented the scheme and obtained 
solutions of /-QCT and of i-QCT for the quantum Fourier transform and a unitary operator that 
does not have an apparent symmetry. The former has a polynomial gate complexity and the latter is 
expected to have exponential one which is based on the fact that a series of generic unitary operators 
has a exponential gate complexity. The time complexity for the former is found to be linear in the 
number of qubits, which is understood naturally by the existence of an upper bound. The time 
complexity for the latter is exponential in the number of qubits. Thus the both targets seem to be 
examples satisfyng the statement above. The typical characteristics of the optimal Hamiltonians 
are symmetry under time-reversal and constancy of one-qubit operation, which are mathematically 
shown to hold in fairly general situations. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Ca, 02.30.Xx, 02.30.Yy 
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I. INTRODUCTION 

Quantum computation is performed by physical pro- 
cesses obeying quantum mechanics. It became one of 
the most exciting field in physics and information the- 
ory after Shor [l[ discovered an algorithm to factorize 
integers which is exponentially faster than any known 
classical ones. In quantum computation theory (QCT), 
as in classical computation theory, the computation time 
is usually defined by the number of elementary steps or 
gates necessary to perform a computation, i.e., to real- 
ize a desired unitary operator. Minimum such number 
is called the gate complexity. We shall call this conven- 
tional QCT as gate-optimal QCT (#-QCT). 

In this paper, we investigate time-optimal quantum 
computation theory (i-QCT) where the computation 
time is defined by physical time. There are two moti- 
vations for £-QCT. The first is to develop (an abstract) 
£-QCT itself which is a physical-time-based alternative to 
g-QCT. Since a quantum computation is a physical pro- 
cess, it is physically more natural to measure the time 
by the physical one. From this viewpoint, the computa- 
tion time in g-QCT can be seen as information-theoretic 
time which is a more abstract or coarse-grained notion 
of time than the physical one. Time optimality is at- 
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tracting growing attention in quantum optimal control 
theory (OCT) mainly in the context of physical applica- 
tions such as control of an atom by an electromagnetic 
field and NMR quantum computation (0, Q and refer- 
ences therein). 

The second motivation, which we stress in the present 
work and state in detail in Sec. Ill B[ is that £-QCT may 
be a useful tool to analyze g-QCT. Finding the gate- 
optimal algorithm is a discrete and combinatorial prob- 
lem, which makes construction of a general theory dif- 
ficult. On the other hand, time-optimal algorithms are 
smooth curves in a certain space which obeys a differen- 
tial equation [U-[6|], typically that for a sub-Riemannian 
geodesic 0] on a manifold. This may allow a general 
theory and approximation methods. Moreover, roughly 
speaking, upper and lower bounds for gate complexity 
can be given in terms of optimal physical time 0, 
Thus i-QCT is useful in the investigation of g-QCT. One 
may become able to calculate gate complexity by calcu- 
lating time complexity. We ask to what extent holds the 
statement that the time complexity is polynomial in the 
number of qubits if and only if so is the gate complexity. 

We will therefore compare the time complexity and the 
gate complexity for some typical examples. We choose 
the quantum Fourier transform (QFT) as an example 
of the target unitary operator of which a fast algorithm 
in the sense of g-QCT (i.e. whose gate complexity is 
polynomial in the number of qubits) is known, while we 
choose a target unitary with no special symmetry because 
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a generic series of target unitary operators is known not 
to have fast algorithms. 

To achieve it, we will make use of an efficient numerical 
method for OCT, so-called Krotov's method. We relate 
f-QCT to fidelity-optimal QCT (/-QCT), and develop 
a Krotov-like scheme for /-QCT by making use of the 
formal similarity of /-QCT and OCT. In the context of 
OCT, similar ideas of replacing time optimality to fidelity 
optimality have been used. 

We will see that both the QFT and the asymmet- 
ric unitary operator satisfy the statement above. Fur- 
thermore, we will find some characteristic behaviour of 
the optimal Hamiltonian, namely, time-reversal symme- 
try and constancy of one-qubit Hamiltonian components. 

In analyzing f-QCT, it is useful to combine numerical 
and mathematical approaches. We will show mathemat- 
ically that the behavior which is found numerically is 
satisfied in fairly general situations. These arguments in 
turn support the soundness of our numerical calculation. 

The organization of the paper is as follows. In Sec. [H] 
we introduce f-QCT as a special case of quantum brachis- 
tochrone and explain our above-mentioned motiva- 
tion more precisely. In Sec. IIII1 we discuss the relation 
between f-QCT and /-QCT. In Sec. II VI we present a 
Krotov-like numerical method for /-QCT. We will show 
our numerical results and extract the properties of the 
solutions of /-QCT and of f-QCT in Sec. El We will 
give a proof of time-reversal invariance and constancy of 
one-qubit components in Sec. (IVI|) . Sec. IVIII is for con- 
lusion. In Appendix [X] we show the monotonicity of the 
numerical scheme of Sec. IIV1 and in Appendix [Bl we give 
a proof of a theorem in Sec. IVII 

We use the units h = 1 . 



II. TIME-OPTIMAL QCT 

In this section, we introduce f-QCT, with some review 
of quantum brachistochrone, and present the motivation 
of the paper in more detail. 



A. Definition 

Let us define the simplest version of time-optimal QCT 
(t-QCT) as a special case of quantum brachistochrone 
for unitary operations [6j, namely, the case in which the 
Hamiltonian H(t) involves only one- and two-qubit in- 
teractions and is subject to a normalization constraint. 
Below is a summary of the formalism Q in this case. 

Quantum brachistochrone for unitary operations is a 
framework to find the optimal Hamiltonian H{t) G T 
which realizes the desired unitary U f up to phase in the 
minimum time T 0, where T is the set of available 
Hamiltonians. Namely, one wants to find the minimum 
T such that there is a unitary operator U(t) satisfying 



the Schrodinger equation 

iU(t) = H(t)U(t), 

and the initial and final conditions, 

U(0) = 1, 
U(T) = e- lx U f , 



(1) 

(2) 
(3) 



where 1 is the identity operator and \ is some real. 

For f-QCT of the system of n qubits, the set T of 
available Hamiltonians consists of self-adjoint operators 



(4) 



(5) 



H{t) = Y,K{t)r a 

a 

with a normalization constraint 



where N := 2" and the basis {r a } consists of cr^/vW 
and <j^/Vn. Here, o^V.'.J™ (1 < ai < • • • < a m sC n and 
ji = x, y, z) denotes the direct product of the Pauli oper- 
ator o~j t on the a;th qubit and identities on the others; for 
example, a^y = cr x ®l® a y ®l® ■ ■ - <£)1. The normaliza- 
tion condition (JS|) can be interpreted physically as finitc- 
ness of available energy in operations, while it is needed 
mathematically for the optimality problem to be well- 
posed Q. The A-dependence of J3]) is for consistency 
under composition of systems [10| ■ The parameter w can 
be interpreted as defining a unit for T. The problem is a 
particular case of linear homogeneous constraints @, Sec. 
III]. 

This is the natural counterpart in f-QCT to the stan- 
dard paradigm in g-QCT where one constructs a desired 
unitary by a sequence of one- and two-qubit operations. 
The parameter lo can be interpreted as defining a unit for 
T. The problem is a particular case of quantum brachis- 
tochrone for linear homogeneous constraints @, Sec. III]. 

In the form of variational principle, f-QCT is to mini- 
mize the action 



S(U, h, V, A) = 



rff 



L s +^-(\h(t)\ 2 -Nu 2 ) 



(6) 



Trt7t(l-P[/)(C7) 



Pu(A) 



where L T ■- y Tr (Huyt (i-Pu)(HU) > 
jj(TrAW)U, L s := 2Re Tr V"t (ijj _ HU) , an overdot 
denotes time derivative. The first term counts the time 
duration, where Lt is unity when the Schrodinger equa- 
tion (fTJ) holds and is invariant under time reparametriza- 
tion t i—)- /(f) [ll|. The second term Ls guarantees that 
the Schrodinger equation (JlJ holds at all times where the 
unitary operator V(t) is the Lagrange multipliers. The 
third term guarantees the normalization constraint ([5]) 
where the real function A(f) is Lagrange multiplier. We 
have adopted an action equivalent to but slightly differ- 
ent in form from that in Ref. Q for better connection 
with the arguments below. 
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We note that the phase of U(i) does not matter in the 
present formulation of i-QCT. In fact, the action (| 15[) is 
invariant under a time-dependent phase change of U(t) 
which can be considered as a gauge transformation [6|. 
Therefore the theory is defined on U(7V)/U(1), and one 
can also think of it as a theory on SU(7V) (by "gauge 
fixing" ) . 

The Euler-Lagrange equations are Q the Schrodinger 
equation (fTJ for U(t), the normalization constraint ([5]), 
the Schrodinger equation for V(t), 



iV(t) = H(t)V(t), 
and the equation determining H(t), 



(7) 



X(t)h a (t) = Tr T a F(t), F(t) := U{t)V\t) + V(t)rf(t), 

(8) 

where X(t) :— X(t)—j^. One must solve these equations 
with the initial and final conditions © and ([3]). 

Let us recall some general features of the system. First, 
F satisfies the simple evolution equation 

iF(t) = [H(t),F(t)], (9) 

Second, the expres- 



which follows from (TTJ, ([7]) and 
sion for X(t) is 



(10) 



which follows from J5]) and ©. Third, this variable X(t) 
is constant in time [12j . 

In the special case when Uf is a one- or two-qubit op- 
eration in \J(N), the solution of i-QCT is given by a Rie- 
nianman geodesic U(t) = e~ lff * on U(JV)/U(1), where 
H is constant 4, 6]. The time T(Uf) is proportional to 
the arc length and depends solely on the eigenvalues e %6 i 
of the relevant SU(4) part of Uf as 



\ 



mm 



E( 



2irmj - x) 2 



(11) 



where \ is a real and m,j are integers. For example, 

the time is given by T([/ CNO t) = T(S) = for 
the controlled NOT gate U C not and SWAP gate S 
whose eigenvalues are (1,1,1,-1), and T 2imax = ^jj- for 
the hardest two-qubit operation whose eigenvalues are 
(1, i, — 1, — i). These can be used for the unit of T(Uf) in 
comparing T(Uf) with G(Uf). 



B. t-QCT, g-QCT and our motivation 

There are some rigorous relations between the gate 
complexity and time complexity. Very roughly speaking, 
one can give upper and lower bounds for gate complexity 
through the time complexity. 



Before introducing the relations, recall that in the sim- 
ple cases, the time complexity T(U) can also be inter- 
preted as the arc length L(U) of the sub-Riemannian 
geodesic connecting 1 and Uf, up to overall multiplica- 
tive constant. The simplest version of i-QCT presented 
in Sec. Ill Al falls into this category, so that j^^ax appear- 
ing below in this paper can also be considered as L ^ 

° -^2, max 

where £2, max is the sub-Riemannian geodesic distance 
between 1 and the furthest two-qubit operation. 
The precise relations are given [9|, esp. Eq. (15)] by 



^2.max ^ 



(12) 



where T(U) is the time complexity to realize the unitary 
17, and G(U, e) is the gate complexity to realize a unitary 
within the distance e from U (measured in the operator 
norm), £2, max is the constant defined in the previous sec- 
tion, and A is some constant. Note that we always have 
G(U,e)^G(U). 

Let us explain the first inequality of (|12|) which is sim- 



T(U) 



pie. For a one- or two-qubit gate U, we have „ 

1 = G(U) by the definition of T2 jlnax . Then, for a general 
U, letting U = U m ■••U\ be the gate-optimal decompo- 
sition, we have Ta Sj=i T(Uj) ^ G(U). However, 
the sum in the LHS, the time cost of this decomposition, 
must be greater than or equal to the time complexity 
T(U) of U. Thus the first inequality of (JT2J) holds. 

The relation (fT2| is suggestive, and one might expect 
that 



G(U)&T(U), 



(13) 



by which we mean that G(U) and T(U) are bounded by 
some polynomial of each other and n. However, it is 
argued that this cannot be true in general @ . Then, one 
may want to ask to what extent (|1 3|) is true in general, 
since (|12[) is derived by a completely general argument. 
In particular, we think that it is interesting to ask what 
are the classes of unitary operators U which satisfy 

1. T(U) w G(U), i.e., T{U) and G{U) are bounded 
by a polynomial of each other and n, 

2. T(U) - G(U), i.e., T(U) is polynomial in n if and 
only if so is G(U). 

The second class contains the first class. Note that 
the first and second conditions above are equiavlent to 
G{U) « G{U,e) and G(U) - G(U,e), respectively, by 
(1121) . The questions above are not easy, but we want 
to develop a basis here which will help answering these 
questions in the future. 

i-QCT itself is a good framework to analyze these ques- 
tions theoretically. However, a drawback is that it is ex- 
tremely hard in practice to derive exact solutions when n 
is large. One way is to appeal to numerical calculations, 
but it is also difficult because one encounters a two-point 
boundary value problem in dimensions rapidly increasing 
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with n. In Ref. 13;], time-optimal solutions were obtained 
up to seven qubits for the like of Hamiltonians appearing 
in NMR quantum computers. That problem involves 2n 
(< 20 for n < 10) functions for the boundary value prob- 
lem. Our t-QCT has 9n(n — l)/2 + 3n functions, where 
the number easily becomes several hundred. To minimize 
the numerical difficulty, we introduce a problem bridging 
t-QCT and optimal control theory (OCT), and make use 
of an efficient numerical scheme for the latter. 



III. FIDELITY-OPTIMAL QCT 

In this section, we relate t-QCT to /-QCT and pre- 
pare for introducing the numerical method in the next 
section. In the context of OCT, similar ideas of relating 
time optimality and fidelity optimality are used (e.g. |3|]). 



B. Relation between solutions to /-QCT and 
t-QCT 

The solution of t-QCT can be obtained from that of 
/-QCT in the limit T -> 1. 

First, /-QCT is equivalent to the problem of minimiz- 
ing the physical time to achieve given fidelity. Suppose 
the solution Uo(t) of /-QCT for given time To, with opti- 
mal fidelity T{U (To), [//) = /, was not the solution of the 
above new problem. There would be U{t) achieving the 
same fidelity in some T < Tq. Then one could construct 
U(t) with F{U{T ), Uf) > / by defining U(t) = U(t) for 
^ t ^ T and appropriately in T < t ^ To. This would 
contradict to the fact that Uo(t) is a solution of /-QCT. 
Thus Uo(t) is the solution of the above problem. The 
converse is shown similarly. 

Second, the above new problem with fixed fidelity / 
yields t-QCT when / — > 1. Therefore the solution of 
/-QCT gives that of t-QCT in the limit T -> 1. 



A. Definition 

Let us define fidelity- optimal QCT (f-QCT) as a 
framework to solve the following problem: given a target 
unitary operator U / and time interval T, find the Hamil- 
tonian H (t) 6 T as a function of time which maximizes 
the fidelity F{U (T), Uf). Here we use the trace fidelity 



T(U(T),U f ) 



(14) 



because we want U(T) to be close to Uf only up to phase. 
The problem is to minimize the action 

S(U, h, V, A) = -NT{U{T), U s f 



dt 



L s + ^-(\h(t)\ 2 -Nu 2 ) . (I.:,) 



The first term is the squared fidelity of the final unitary 
U(T) with respect to the target Uf. The second and 
third terms guarantee the Schrodinger equation (fTJ) and 
the normalization constraint (JSJ, respectively, as in the 
case of t-QCT. The action (|T5|) is invariant under the 
(time-dependent) change of the phase of U (t) so that the 
theory is defined on U(AT)/U(1). 

The Euler-Lagrange equations yield jT|), (0, flEJ, and 



V(T) = -U f TrU}U(T). 



(16) 



Therefore all the equations (fl"1)- (fTU|) hold except that §S§ 
is replaced by ([To]) . As in the case of t-QCT, A(t) is 
constant in time. 

The solution gives the maximal fidelity F(U(T),Uf) 
for given T. /-QCT is important on its own when one 
discusses the tradeoff between speed and error of compu- 
tation, here we use /-QCT to bridge t-QCT and OCT. 



C. Formal similarity of /-QCT and OCT 

Let us see the formal similarity between /-QCT and 
OCT. 

Our action (| L5[) defines an exact fidelity optimality 
problem with constraints. If we replace A(t) with a given 
constant or a given function, the term j n t ne in- 

tegral can be interpreted as a penalty term while the con- 
stant term — h^lM^i m ay be dropped. Then the problem 
is to minimize a combination of the fidelity and penalty 
terms with their weights specified by A(t), which is a typ- 
ical problem in OCT. This is the formal relation between 
/-QCT and OCT. 



IV. A KROTOV-LIKE SCHEME 

In this section, we shall define an efficient numerical 
scheme for /-QCT by modifying Krotov's method [13j in 
OCT, making use of the similarity between /-QCT and 
OCT in Sec. HHCl 

In what follows, the functions h a (t) and A(t), respec- 
tively, are h a (t) and A(t) calculated in the middle of an 
iteration cycle of the scheme, and Eqs. (O, and ©, 
respectively, denote Eqs. (0, and (|SJ) with h a (t) (in- 
cluding that in H{t)) and A(t) being replaced by h a (t) 
and A(t). The scheme is as follows. 

(i) Prepare a seed Hamiltonian components h a (t), ^ 
t < T. 

(ii) Set U(0) = 1 and evolve U(t) from t = to t = T 
by (fTJ) with the Hamiltonian H (t) = ^ a h a {t)r a . 

(iii) Set V{T) by (16| and evolve (V(t),h a (t),X(ff) 
backward in time from t = T to t = by Eqs. j7|), 
© and (flU)) . while U(t) is treated as a given function 
and is not updated. 
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(iv) Set U(0) = 1 and evolve (U(t), h(t), X(t)) forward 
in time from t = to t = T by Eqs. (JTJ, © and (|TU|) . 
while V(t) is treated as a given function. 

(v) Repeat steps (iii) and (iv) until the variables con- 
verge; the final h a (t) defines the optimal Hamiltonian 
and the final T(U(T),Uf) gives the maximal achievable 
fidelity in T. 

Note that the multiplier X(t) converges to an unknown 
constant only after the convergence, which is in contrast 
to original Krotov's scheme where A is a given constant 
parameter of the problem. The constancy of A(t) can be 
used for a convergence check of the present scheme. The 
present scheme is also different from Krotov's method for 
the case of a fixed reference energy (see, e.g. [Hj]). In 
the former, the variable X(t) is a Lagrange multiplier and 
the normalization (JS|) must be satisfied at all times, and 
X(t) is updated in the scheme. In the latter, X(t) is a 
given function which determines the amount of penalty 
imposed on the error in ([S]), and X(t) is not updated in 
the scheme. See also discussions in Sec. IIIIBl 

An important property of our Krotov-like scheme is 
monotonicity, which is necessary for the scheme to be 
useful. We give the proof in Appendix [Al 



RESULTS 



unitary operator U / with a generic image of a fixed state 
vector, Uf\0), has the gate complexity G(Uf) ~ 2™ [l9j . 

(n) 

Our concrete choice is I7asym which is, in the matrix form, 



w asym 



7o a 



7ia 7 2 a 
7i/3i 7 2 ai 







7o a* 



72/#2 






7n-l a l 
7n-l«2 



7n-l/3n- 



(19) 



where a k := (k + l) 1/3 e iv/ ^, and fa € C and jj > 
are determined by orthogonalization and normalization, 
respectively, of the columns. The state vector U^y m \Q) 
is given by the first column of the right hand side of 
(|19[) which does not have apparent symmetry. One can 
therefore expect that the gate complexity is exponential 
in n, G{UjS$m) > 2". 

For each given time T, the convergence of the scheme 
is checked by the convergence of the fidelity, and the con- 
stancy of A(i) explained in Sec. IIVI 



We implemented the numerical scheme presented in 
Sec. IIVI and performed calculations for /-QCT. We chose 
two examples for the target unitary operator Uf. 

The first is the quantum Fourier transform (QFT) 

C/qft> defined by 

k=N-x 

U<$ I \x) = -±= ^ lkx,N \k), x = 0,l,..,N-l, 

fc=0 

(17) 

which has the gate complexity G(Uf) polynomial in n. 
In fact, a simple efficient algorithm (e.g. is given by 

^QFT = ^nzjti— na ' ' ' ^lnUi ■ ■ ■ U ni 

Uj '■= Rj,n-j+l,n ■ ' ■ R2.n-j+l,n-j+2W n -j + i, (18) 

where Wj is the Walsch-Hadamard gate W :— 



1 

V2 



1 1 

1 -1 



applied on the jth qubit, Rj t km is the 2tt/2 j - 



phase shift gate Rj := 



on the fcth qubit con- 



"1 

e 2 "/ 2J 

trolled by the mth qubit, Sjk is the SWAP gate S on 
the jth and fcth qubits, and n 2 := [n/2\ (integer part of 
ri/2). The number of gates of this construction is n(n + 

l)/2+[n/2j so that G(U$ T ) < n(n+l)/2+ [n/2\ n 2 . 

The second example of the target Uf is chosen so that 
we can expect that Uf has the gate complexity G(Uf) 
exponential in n. To do so, we pick a Uf which does not 
have any apparent symmetry because a generic unitary 
operator Uf has the gate complexity G(Uf) « 4" and a 



A. Fidelity-time relation 

Let us show the fidelity-time relation, namely, the max- 
imal achievable fidelity F(U(T), U f) in time T. For sim- 
plicity, we only show the case of the QFT (Fig. [T]). 

Though our scheme in Sec. IIVI is monotonic, there is a 
possibility of the calculation is trapped by a local min- 
imum of the action S. We did the following two things 
in order to find the global maximum J- . One is simply 
preparing many random seeds h a (t) for each T. Another 
is making use of the continuity of the solutions. Namely, 
we prepared many random seeds h a (t) for some fixed T. 
Then we used the solution for T as the seed for a nearby 
T, and find continuous branches of locally optimal so- 
lutions. This "output recycling" turned out to be often 
more powerful than merely increasing the number of ran- 
dom seeds for every T . We observed crossovers of those 
branches. Only the branches of the largest T contribute 
to the curves in Fig. [TJ 

We observe the following, which may be characteristic 
of the QFT. 

(a) The odd and even qubits seem to make a pair (4-5, 
6-7, and 8-9) for n > 4. 

In other words, the time-optimal solutions seems to 
split into the series of odd n and that of even n, which 
may be useful in the future mathematical analysis of the 
time-optimal solutions of the QFT. 
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FIG. 1: (Color online) (Upper) Optimal achievable fidelity 
in given physical time for the n-qubit QFT in /-QCT. The 
numbers beside the curves denote n. The n = 10 result is pre- 
liminary. One can see property (a) in Sec. IV Al (Lower) The 
same graph with the vertical axis being 1 — T in logarithmic 
scale. The points and the numbers on the lower axis indicate 
the time complexity T(Uf) estimated by the fidelity — > 1 limit 
for each n. 



B. The limit T -> 1 
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FIG. 2: (Color online) Estimation of T(Uf) from the limit 
J" ->■ 1. In the n = 5 QFT example, y := 1 - T{U(T), U f ) 
is fitted by i/ = 0.735(1.81 - x) 2 ' 84 , where a; := T/T 2 ,max, 
giving an estimate T((7/)/T2,max — 1.81. We used the data 
0.002 sC y 0.01. 



C. Time complexity as function of n 



Let us discuss the behaviour of the time complexity as 
a function of number of qubits, n. 



Fig.[3]shows the time complexity T(C/qf T ) in i-QCT as 

a function of n, where T(C/qf T ) is obtained by the value 
of T on each curve in Fig.[JJin the limit J- — > 1. For n = 1 



and n = 2, one has the exact values 



T 2 - 



-I 

2 I 



4 



2 



0.8944 and 



8/4 



vTi 

2^/5 



0.7416, 



^2 , max 

calculated from dill) , with which our numerical results 
agree well. 



We estimate T(U) from the limit F — > 1 of the solu- 
tions to /-QCT in Sec. IVH 

In the limit, we have two sources of error. One 
is a natural numerical error which makes the fidelity 
J-(U(T),Uf) saturate below unity Another is that if 
T is larger or very close to the time complexity T(Uf), 
the solution of /-QCT for the physical time T begin to 
"take a roundabout route" before reaching U(T). 

With this behaviour in mind, we estimate the time 
complexity by a nonlinear fitting of the fidelity-time 
curve around ^"({/(T), Uf) <; 0.99. Fig. [2] shows an ex- 
ample of the estimation, in the case of the n = 5 QFT. 
We fit y := l-T{U{T), Uf) by y = a(b-x) c , where x := 
T/TVmax- The estimation is given by T(J7/)/r 2 ,max — b. 
The error in the time complexity T(j7)/T2 )max is about 
0.1 in the case of QFT, and it is about 0.1 or 0.2 in 
the case of asymmetric unitary operator. This does not 
change the conclusion of the subsequent sections. 



We observe the following property from Fig. |3l 

(b) The optimal time T([/qp? T ) is linear in the number 
of qubits, n. 

The line in Fig. [3] is the result of a linear fitting, which 
is T(Uf)/T 2 , max = 0.32n + 0.27. We used the data n > 2 
in the fitting because the behaviour of n $J 2 and that of 
n ^ 2 should differ due to the nature of T allowing only 
interactions involving two qubits or less. 

Property (b) is in good contrast to the number of 
gates, 0(n 2 ), of the known efficient algorithm (ITS)) for 
the n-qubit QFT. However, it can be understood natu- 
rally. Since 



T(W) _ 



_2_ _ 

T2,max v5' ^2, max y 5 T2,max 

3 from (|11[) . we have the physical time cost T'( n > 



T(S) = 
T*2 max 



and 
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of the construction (|18|) is [Vi 

rp/(n) 



-^2, max 
1 



To 



nT{W) + J2( n - J + l ) T (Ri,km) 

3=2 



T(S) 




(20) 
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which is bounded (from above and below) by a linear 
function of n. The time complexity T(Uqp T ) is several 
times smaller than T'( n > (except for n = 1 when they 
coincide). The significance of J 1 ''™) is that it is a rigorous 
upper bound for the time complesity, T(Uq^ t ) ^ T'( n ) ; 

which implies that T(Uq^ t ) is at most linear in n. This 
strongly supports property (b) and the correctness of the 
numerical calculation. 

Fig. 0] shows the time complexity T(Uf) as a function 
of n for the case of the asymmetric target (7i™y m . 

We observe that 

(b') The optimal time T(llisym) is exponential in the 
number of qubits, n. 

As in the case of QFT, we used the data n > 2 to 
fit by a function. The time complexity T(l7asy m ) is well 
fitted by an exponential function as T(Ui2ym) /l2 jma x = 
0.20 x 2 S2n . It is suggested from the numerical result 
that U$m is in the class T(U) ~ G(U). 

To conclude, it is suggested that both the QFT and 
the asymmetric unitary are in the class T(U) ~ G(U), 
which is polynomial in n for the former and exponential 
in n for the latter. [Note that in the polynomial case, 
T{U) ~ G(U) implies T(U) « G(U)]. 



FIG. 3: (Color online) Optimal time T{Uj) for realization of 
the QFT as a function of the number of qubits, n, with a linear 
least squares fitting for n ^ 2: T(U f )/T 2 ,ma,x = 0.32n + 0.27. 
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FIG. 4: (Color online) Optimal time T for realization of 
the asymmetric unitary Ua"-ym as a function of the number 
of qubits, n, with a linear least squares fitting for n ^ 2: 
T(Uf)/T 2 ,ma,x = 0.20x2°- 82n . Note that, in contrast to FigEJ 
the vertical axis is in logarithmic scale. 



D. Behavior of the time-optimal Hamitonian 

Let us analyze the behavior of the optimal Hamiltonian 
H(t) in f-QCT. 

We shall say that an element of su(iV) is symmetric 
(or antisymmetric) if it is so in the standard matrix rep- 
resentation. In particular, r a is symmetric (antisymmet- 
ric) if it contains even (odd) number of a y ; for example, 
ayy/y/N is symmetric and a^/y/N is antisymmetric. 

Fig. [5] is the behavior of the Hamiltonian H(t) for the 
4-qubit QFT with T ~ I, which can be considered as 
the solution of i-QCT. The components h a are catego- 
rized into the four according to: whether r a is one- or 
two-qubit interaction, and whether r a is symmetric or 
antisymmetric. 

The results suggest the following. 



(c) The components for symmetric (antisymmetric) r a 
is symmetric (antisymmetric) under time-reversal t i— ► 
T — t, and 

(d) one-qubit interaction components are constant. 
Note that (c) and (d) imply that one-qubit antisymmet- 
ric components vanish. 

The same properties are satisfied by the n — 5 QFT. 
For the n = 3 QFT, (c) does not hold but (d) does. 

(n) 

For the asymmetric target [/asym, we do not observe 
property (c), the time reversal invariance. However, we 
do obverve property (d), constancy of one-qubit compo- 
nents, also in these cases. Fig.|5]shows the n — 3 example 
of the behavior of the optimal Hamitonian. 

These properties will be discussed from a theoretical 
point of view in Sec. IVIl 
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FIG. 5: (Color online) Behavior of the optimal Hamilto- 
nian for the 4-qubit QFT with T/T 2 ,ma X = 1.48, T = .9967. 
Shown are the components with respect to the generators 
which are (Is) one-qubit symmetric, (la) one-qubit antisym- 
metric, (2s) two-qubit symmetric, and (2a) two-qubit anti- 
symmetric, in the sense of Sec. IV CI The solution is time- 
reversal invariant. The one-qubit components are constant in 
time. 
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FIG. 6: (Color online) Behavior of the optimal Hamiltonian 
for the 3-qubit asymmetric unitary U^ym with T/T2, m ax = 
1.03, J- — .9997. Shown are the components with respect to 
the generators which are (Is) one-qubit symmetric, (la) one- 
qubit antisymmetric, (2s) two-qubit symmetric, and (2a) two- 
qubit antisymmetric, in the sense of Sec. IV CI The solution 
is not time-reversal invariant. The one-qubit components are 
constant in time. 



VI. MATHEMATICAL JUSTIFICATION OF 
THE BEHAVIOUR OF THE TIME-OPTIMAL 
HAMILTONIAN 

The temporal behaviour of the optimal Hamiltonian 
in £-QCT found in Sec. [V] property (c) for some of the 
QFT and property (d) for the QFT and the asymmetric 
unitary, are not peculiar to the case of those target uni- 
tary operators. They in fact can be proven under a fairly 
general condition. These results also serve as evidences 
of reliability of the numerical calculation. 

A. Time reversal invariance 

Let us show (c), the time- reversal symmetry found in 
some of the time-optimal solutions in Sec. [Vj 

Let us define the time reversal 
(U Tev (t), iT rev (t), Vrev(i), Arev(i)) of the set of vari- 
ables, (U(t),H(t),V(t),X(t)), by 

U ICV (t) :=U*(T-t)U T (T), 

H Tev (t) :=H*(T-t)=j: a h a (T-t)T* a , 

Ke V (i) :=V*(T~t)U T (T), 

Arcv(i) := A(T — t), (21) 

where the superscript asterisk denotes the complex con- 
jugate and the superscript T denotes the transpose. We 
have the following, whose proof is given in Appendix [B] 

Theorem 1. Let the target Uf be symmetric 
up to phase. If (U(t), H(t),V(t), \{t)) is a so- 
lution of t-QCT for Uf with optimal time T, 



so is (U iev (t),H lev (t),V iev (t),X le v(t)). In par- 
ticular, if the solution is unique, it is invari- 
ant under time-reversal, (U(t), H(t),V(t), A(t)) = 
( U ICV (t), H Iev (t),V Tev (t), A rcv (t ) ) . 

The QFT is a symmetric target. The theorem justi- 
fies the observed symmetry property (c) of the QFT for 
n =/= 3. Convergence of ramdomly chosen initial Hamilto- 
nians to a single time-symmetric solution suggests that 
the optimal solution is unique in those cases. The 3- 
qubit case did not show the symmetry, and consistently 
we observed many optimal solutions H{t). 

The asymmetric target unitary is not a symmetric op- 
erator and the numerical solutions did not show time- 
reversal invariance. 

We remark that Theorem [T] holds not only in the 
present version of t-QCT but also in any quantum 
brachistochrone with T* = T. 



B. Constancy of one-qubit components 

Let us prove (d), which turns out to hold in general. 

Theorem 2. The one-qubit part of the Hamiltonian for 
any solution of t-QCT is constant in time. 

Proof. Let Qj be the space of j-qubit operations. They 
satisfy the following commutation relations [y, Sec. V]: 

[Sj,0k] = 0\j-k\+l © 0|i-fc|+3 © ' ' ' © 9j+k-l, (22) 
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where j,k ^ 1 and gj := for j > n. From ©, we can 
decompose F as 

F = XH + F' (23) 

with H = J2j=i 2 Hj an{ i ^' = Sj>3 ^j'i where an in- 
teger subscript j denotes the projection to Qj. Then 
it follows from (j2"2"]l and A = const, that the equation 
iF = [H, F] can be written as 

iHi = 0, 
iXH 2 = [H 2 ,F 3 } 2 , 
iF 3 = [Hi,F 3 ] + [H 2 ,F A ] 3 , 

(24) 

This implies H± = const, for any target unitary Uf. □ 

Since constancy of one-qubit components is a quite 
general feature, it is a useful criterion of the convergence 
of the numerical scheme. 



VII. CONCLUSION 

We investigated the simplest version of t-QCT, where 
the computation time is defined by the physical one and 
the Hamiltonian contains only one- and two-qubit inter- 
actions. This version of i-QCT is also considered as op- 
timally by sub-Riemannian geodesic length. 

Motivated by the relations between time complexity 
and gate complexity (IT21 . we aimed to pursue the possi- 
bility of using time complexity as a tool to estimate gate 
complexity, and asked the following question: to what 
extent is true the statement that time complexity is poly- 
nomial in the number of qubits if and only if so is gate 
complexity. In particular, we want to identify the classes 
of unitary operators U that satisfy T(U) ~ G(U) and 
T(U) ~ G(U), by which we meant T(U) and G(U) are 
bounded by polynomial of each other and n, and T(U) 
is polynomial in n if and only if so is G(U), respectively. 

For this program, we introduced an efficient Krotov- 
like numerical scheme by making use of the relation be- 
tween i-QCT and /-QCT and the formal similarity of the 
latter to OCT, and showed its monotonic convergence 
property. 

We chose the quantum Fourier transform as an ex- 
ample of the target with polynomial G(U) and a uni- 
tary operator without symmetry that is expected to have 
exponential gate complexity. We obtained the fidelity- 
time relation, time copmlexity T(U), The time complex- 
ity of the QFT is found to be linear in the number of 
qubits. The time complexity of the target without sym- 
metry is exponential in n. These results suggest that the 
QFT and the asymmetric target are both in the class 
T(U) « G{U), and that T(U) ps G(U) is linear in n for 
the QFT and is exponential in n for C" _1 -NOT, respec- 
tively. This supports the usefulness of time complexity 
as a tool to estimate gate complexity. It is also suggested 



that a polynomial-gate algorithm does not exist indeed 
for the asymmetric target llisym, because T(U)/T 2jmaDl is 
the absolute lower bound of G(U). 

We also found two characteristics of the optimal Hamil- 
tonian H(t) of f-QCT. One is symmetry under time re- 
versal and the other is constancy of one-qubit operation, 
which are mathematically shown to hold in fairly general 
situations. 

A natural extension of this work is to push forward 
with the program above by comparing the time complex- 
ity (or equivalently, the arc- length of the sub-Riemannian 
geodesies) and the gate complexity for other unitary op- 
erators. Another direction is to consider other variants of 
i-QCT. An example is t-QCT which allows only nearest 
neighbor interactions in a lattice, and another is t-QCT 
where only the time spent in two-qubit interactions is 
counted and that spent in one-qubit interactions is ne- 
glected. 

We would like to stress that interplay of numerical and 
mathematical methods, of which an example was the ar- 
gument given in Sec. I VII is important in analyzing t- 
QCT. We hope that our method will lead to new under- 
standing about the power of quantum computation. 
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Appendix A: Monotonicity of our Krotov-like 
scheme 

Let us show the monotonicity of the Krotov-like 
scheme given in Sec. IIVI 

Let AS be the change of S in the cycle (iii)-(iv). 
We would like to show AS ^ 0. Let the variables 
(U(t),h(t),V(t),X(t)) be 

(U hl (t),h in (t),V ia (t),X in (t)) 

after step (iv) which will be the inputs to a new cycle 
(iii)-(iv). They satisfy (JTJ, ©, ©, ©, but not © or 
([To). Let the variables be 

{U in (t),h(t),V out {t),X(t)) 

with V out = Vi n + SV and h — h m + 8h after step (hi). 
They satisfy ©, ©, ©, (JE}, but not © or ©. Let 
the variables be 

(U out {t),h out (t),V out (t),X out (t)) 

with [/ ou t = Ui n + SU and h out — h + Sh after step (iv). 
They satisfy ©, ©,©,©, but not © or JIB]). 
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The change AS in the action (fT5|) after one cycle (iii)- 
(iv) is given by 



NAS = - N 2 F{U out (T), Uff + N 2 F(U m {T), U f ) 2 

ljSU{T) - \TrU} 



2ReTr u} n {T)U } TrU)5U{T) - | Tr u)8U{T)\ 2 , 

(Al) 



because both (£/j n , h in ) and (U on t, h out ) satisfy ([T]) and (O 
and make the integrand in (I15[) vanish. The first term on 
the RHS of ED is 



2Re Tr iV^ ut (T)SU(T) = - I dt 2Re Tr i(V* ut 6 U)' 

Jo 

pT _ 

dt(X out h out ■ 5h + Xh- 8h) 



1 



dt(X out \Sh\ 2 + X\8h\ 2 ) ^ 0, 



(A2) 



where the dot denotes the inner product and we have 
used 



V out (T) = -U f TrU}U in (T), 



5*7(0) = 0. 
i5U = SHU out 

iV out = HV ont , 



- SHU in + HSU, 



KutK ut = 2ReTrr Q [7 out t/ t ut , 



(A3) 



which follow from the conditions satisfied by each of the 
variables, given in the previous paragraph. One can see 
the last equality in (|A2[) by noticing that (h out + h)/2 is 
orthogonal to Sh because \h out \ = \h\, and so forth. 

We conclude that our scheme is monotonic because 
ED and G2J) imply AS ^ 0. Note that SU etc. above 
have not been assumed to be small. 



Appendix B: Proof of Theorem [T] 

It is easily verified that if (U (t) , H (t) , V (t)) is a solu- 
tion of CD)-© and so is (U iev (t), H mv (t), V iev (t)). 
Eq. ([1]) is seen by 

iif rcv (t) = -iU*(T - t)U T {T) 

= H*(T- t)U*(T - t)U T (T) = H iev (t)U iev {t). (Bl) 

Eq. (SJ) follows from U lev (Q) = U*(T)U T (T) = 1; Eq. 
© follows from 

?7 rcv (T) = U*{0)U T {T) = e- ix Uj = eT lx XJ s . (B2) 

Eqs. @j and © follow from r* = ±T a . Eq. (|8} is 
equivalent to 

A(t)H(t) = J! r Q Tr r a (U(t)V^ (t) + V{t)rf (*)) , (B3) 



The time reversal satisfies the same equation (|B3j) be- 
cause 

X Iev (t)H Tev (t) = A(T - <)/7*(T - t) 

= J] r a Tr r Q (C/(T - t)^ (T - i) + V^(T - *)C/ + (T - t)) 

a 

= T a Tr r a (C/ rcv (t) V; f ov (t) + V rcv (t)uL (t)) , (B4) 



where we have used r* = ±r a (the signature depends on 
a). Therefore the time reversal of the solution of i-QCT 
with the target Uf and the time T is a solution of the 
same problem. 
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